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Abstract 

Some time ago Kaplan proposed a new model for the description of chiral fermions on 
the lattice by adding an extra dimension for the fermions. A variant of this proposal was 
introduced by Shamir and can be used to describe vector-like theories in even dimensions. 
We used this model for the simulation of the massive Schwinger model at different gauge 
couplings. The prediction that the fermion mass gets only multiplicative renormalization 
was tested and verified. 



^ 1 Introduction 

t > : 

Some years ago Kaplan suggested a new method simulating chiral fermions on the lattice |]]. 
He started with a vector-like odd- dimensional theory and a mass term that depends on one 
space coordinate and has the form of a domain wall. He found a zero-mode bound to this 
domain wall. From the point of view of the even-dimensional domain wall, this zero-mode is 
a chiral fermion. The zero-mode exists up to a critical momentum, so that the chirality is 
a low energy effect. Anomalies in gauge currents arise from the flow of charge into the extra 
dimension. In the case of a finite extra dimension which has the topology of a circle and periodic 
boundary conditions an additional anti-domain wall with a chiral fermion of opposite handness 
is generated. 

A similar but technicaly simpler approach was made by Shamir [[J. Instead of an infinite 
extra dimension and a mass that is a function of this dimension, he used a semi-infinite coor- 
dinate with free boundary conditions and a constant odd- dimensional mass. In the case of a 
finite lattice this corresponds to cutting the links at both domains in the original model and 
dropping the side with a negative mass. Therefore, in numerical simulations one has to take 
only half as many points to obtain similar results. The zero-states with different chirality are 
located at the boundaries. A connection of these two boundaries by a link of strength m will 
form a Dirac fermion from the two Weyl fermions whose mass is proportional to m, if the 
exponentially small overlap between the surface states can be neglected. This model can be 
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used to describe vector-like theories in even dimensions. Unlike the case of chiral gauge theories 
H in the vector-like model there is no problem to couple (2n + l)-dimensional fermions to a 
gauge field in 2n dimensions. It was shown that perturbative corrections to the quark mass are 
proportional to m, so that it undergoes only multiplicative renormalization to all orders ||. 
Therefore, there is no fine tuning problem that appears for Wilson fermions. This model has 
in the massless case a very mild breaking of the axial symmetries. Moreover, the correct chiral 
limit can be proven M. 

We examined the case of (2 + l)-dimensional fermions coupled to a two-dimensional U(l) 
gauge field, i.e. the massive Schwinger model 0, §[. To test the prediction that the fermion 
mass gets only multiplicative renormalization we studied the behaviour of the pseudoscalar 
isotriplet ('pion') mass if m tends to zero. 



2 Free boundary fermions 



We start our discussion with the free fermion case in 2 + 1 dimensions and a lattice of size 
Q = L 2 ■ L s . The action is given by 
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where (f) = X^ =1 cr^d^ and A = Y7u=x °a« are ^ ne P aim matrices, r denotes the Wilson 
parameter and Pr_l = (l±er 3 )/2 are the projection operators. A lattice point is given by (x, s), 
where s labels the extra coordinate. In the following pages you will note that the lattice spacing 
is taken to be a = 1. The lattice derivation and Laplacian A^ for \x = 1,2 are defined as 
usual and 



, , ^x,2 S = 1 

2 { ®x,s+i - *x,«-i 2 < s < L s 

~ S = L S 



(2) 



^x,2 ~ 2^,1 



-x,L a -l 



S = 1 

2 < s < L s - 1 
s = L s 



(3) 



i.e. the case m = corresponds to open boundaries in the third direction. For r = 1 and 
m — 1 the model corresponds to the topology of a circle and periodic boundary condition in 
s, for r = 1 and m = — 1 it describes the antiperiodic case. Thus, in these cases the (2+1)- 
dimensional boundary fermion model is equivalent to a three-dimensional model with Wilson 
fermions and the corresponding boundary conditions. 
The action ([!]) leeds to the Hamiltonian 



H = o 2 



iaxh + a 3 d 3 + M+ 1 - (-kl + A 3 ) + m {S sA 5 s ^ Lb P l + 8 SjLb 8 sI>1 P r ) 



(4) 



k = sm(k) , k = 2sin(fc/2) . 

Its eigenvalues occur in pairs corresponding to particle and anti-particle and were already 
examined for m = In this case the chiral zero- modes are confined to the two boundaries. 



2 



L=10 M = 0.9 r=l 

s 




Figure 1: Part of the energy spectrum for a lattice extension of L s = 10 and a 
(2 + l)-dimensional mass M = 0.9. The three different types of lines characterize 
the different values of the mass parameter m. 



The wave functions have the form of plane waves in the two-dimensional space and decay 
exponentially with (1 — M) in s. The chiral zero-modes exist for momenta k\ below some 
critical momentum k c . For < M/r < 2 there is one zero-mode at each boundary. For 
m ^ 0, an infinite lattice in x\ direction and a finite third direction one can calculate the 
lowest eigenvalues in the limit where m 2 and k 2 are small, i.e. m 2 , k 2 <C M 2 . If one choose L s 
large enough, so that the exponentially small mixing of the two chiral modes can be neglected 
((1 - M) Ls < m), the result is 

E 2 = h 2 + m 2 M 2 {2-M) 2 . (5) 

Therefore, the mass of the light Dirac fermion is 

fn = mM(2 - M) . (6) 

In Fig. [I] we sketch the dependence of the energy spectrum for the two lightest states on the 
mass. One can observe that the lowest eigenvalue has indeed the predicted form. 

3 Interacting theory 

The interacting theory is defined by coupling the fermions to a two-dimensional gauge field, 
i.e. U^x^)^ = U XfJi for /x = 1, 2 and U(^ x ^ s )3 = 1. This leads for r = 1 to the following gauged 
fermionic part of the action 

S F = (*,Q[E/]tf) 
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= E EU E ( X + ^) + (! " ^m**,.] + (M - 3) 

X S=l I ^ /i=l J 

+ El E [* a! ,.^i*x I -+i + *x,.+iPl**,.] + m [* x ,l.-Pr*x,i + *x,l P L#x,iJ ) , (7) 

x I s=l J 

where £7^ denotes the U(l) gauge variable. Apart from the unconventional sign of the mass 
term, the first term is the usual two-dimensional action. The preferred boundary conditions 
are periodic in space direction and antiperiodic in time direction. The gauge part is given by 

S G = /3£(1-Re?7p) , (8) 
p 

where Up is the product of four C/'s around a plaquette and j3 = 1/g 2 . The whole action is 
given by the sum of the two parts: 

S = S F + S G . (9) 

The case m = should yield to the massless Schwinger model by taking first the limit L s — > oo 
and then /3 — > oo. 

The definitions of the vector and axial currents as well as their divergence equations can be 
found in || ||. Here we just give a definition of the axial transformation (for even L s ): 

S\^ x ,s = iq(s)aiip X;S , (10) 

S A^x,s = -^Q(s)^ x>a <Ti, (11) 

where 

/ n _ / 1 1 < s <L s /2 
q{3) = { -1 LJ2 <s<L s ■ ^ 

The non-invariance of the action under this transformation results from the coupling between 
the layers s = L s /2 and s = L s /2 + 1, respectively, s = 1 and s = L s for m ^ 0. 

For the numerical simulations we used the Hybrid Monte Carlo algorithm with pseudo- 
fermions fL0 |. This requires a positive definite determinant and therefore the flavour duplication 
of the fermion spectrum. The fermion matrix for the first species ist Q[Z7], for the second Q [£/]'. 
If we define the reflection 1Z by 

R = M' ,,,. . : s , (13) 

the fermion matrix satisfies 

a 3 nQ[U]a 3 n = Q[U]^ . (14) 

This yields to 

det Q[U] = detQ[f/] t . (15) 

Thus, the simulated model describes the two flavour massive Schwinger model. 
The topological charge is given by 

Qtop = ^E^> Q top e^, (16) 

^ X 

exp (iF x ) = U xA U x+l2 U* x+U U* x>2 , F x E [-vr, vr) . (17) 

The different topological sectors are separated by a potential barrier of height 2(3, so that the 
tunneling is suppressed. We restricted the simulations to the topological trivial case (Qtop = 

1 Contributions of the non-trivial sectors were examined in jllj] 
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Figure 2: Correlation function for the pion for (s, s') = (1,L S ). The error bars 
for m > 0.1 were omitted because, they are too small for a visualization (smaller 
than 3% of the measured values). The curves shown are the best fit in the interval 
13 < t < 16. 



i.e. we start with U Xfl = 1 and controlled that the topological charge doesn't change. The lightest 
particles in this model are a pseudoscalar isotriplet with mass Mi- satisfying 

Mi- 




>c hf > ( 18 ) 




2.07 . (19) 

This is the analogue to the pion in QCD and becomes massless as the fermion mass m goes to 
zero. The matrix inversions were performed by the conjugate gradient algorithm [|l^ . 

It can be shown that in the limit L s — ► oo corrections to the fermion mass m are proportional 
m, so that the fermion mass in the boundary model gets only multiplicative renormalization. 
The reason is that in weak coupling perturbation theory, corrections to the inverse fermion 
propagator are exponentially suppressed with \s — s'\. Therefore, the zero-modes can't mix, 
except for an exponentially small factor, and a mass term can't develop. Thus, the chiral limit is 
attained if m goes to zero. We used the massive Schwinger model with two degenerate fermion 
species to test non-pertubatively the prediction that the fermion mass tends to zero for m — > 0. 

The simulations were performed at M — 0.9, so that at tree level one has m = 0.99 m. 
For the gauge coupling the two values (3 = 0.25 and (3 = 8.0 were examined. In the last case 
the renormalization of M is small. It follows that the exponentially decay in the additional 
direction is fast and the exponentially contribution to the mass are even for small values of L s 
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Figure 3: Pion mass as a function of m at (3 — 8.0. 



strongly suppressed. Therefore, only small extensions of L s are needed. As a check we used 
lattices of size Q = L x ■ L t ■ L s = 12 • 24 ■ 10 and Q = 32 2 ■ 4. For the mass parameter m we have 
chosen the values m = 0.0, 0.1, . . . , 0.5 in the first case and m = 0.01, 0.05, 0.10, 0.15, 0.20 in 
the second case. At (3 = 0.25 we examined m = 0.0, 0.1, 0.2, 0.3 for a 12 ■ 24 • 10 lattice. 
The pion mass Mi- is calculated from the correlation function 

A 7r s ,(*) = (<,(0)7r s ° s ,(t)) , (20) 



(21) 



where (/, J 3 ) p = (1,0) and x = (x, t). For (s,s') 



and 



(L s , 1) this yields to 



E 



(0,0),(x,t);l,L s 



Q[U]- S 



(0,0),(x,t);L s ,l 



(22) 



For every lattice size we performed 10000 sweeps. Every 100th sweep we calculated QlU]^ 1 , the 
correlation function and the pion mass separately on every configuration. The statistical error 
was estimated by binning. The autocorrelation times for the pion mass were of order (9(100) 
sweeps. The maximum value of approximate 500 sweeps was obtained for (3 = 8.0, Q = 32 2 • 4 
and m = 0.01. 

In fig. the different correlation functions for the 32 2 • 4 lattice and (s, s') = (1, L s ) are 
plotted. From this we extracted the pion mass by a least square fit by a single cosh. The 
values are shown in fig. |3| together with the masses of the other lattice. One can see that the 
values of the two masses for m = 0.20 coincide within the statistical errors. For small values 
of m finite size effects appear. An estimate for the value of the critical mass by polynomial 
interpolation gives \m c \ < 0.01, i.e. no deviation of a multiplicative renormalization for the 
mass is observable. Since we are in the weak coupling region renormalization effects are small 
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Figure 4: Pion mass as a function of m at (3 = 0.25. The dashed line indicates an 
ansatz proportional m 2 / 3 . 



anyway. The renormalization of the fermion mass can be obtained by comparing the data with 
the results of other simulations and perturbation theory. The result is fn = 0.90(5) m. 

Renormalization effects at (3 = 0.25 are larger. Therefore the simulations at this point are a 
better check for the prediction that the fermion mass gets only multiplicative renormalization. 
As before there is no hint for a nonzero critical mass (fig. §). Moreover m/g is small enough 
for a comparison with perturbation theory. The results are in agreement with theory and lead 
torn = 0.60(2) m. 

Acknowledgments 

I thank Istvan Montvay and Karl Jansen for helpful discussions and comments. Especially I 
benefit from Elfie's lectures. 

References 

[1] D.B. Kaplan, Phys. Lett. B288 (1992) 342. 
[2] Y. Shamir, Nucl. Phys. B406 (1993) 90. 

[3] M.F.L. Golterman, K. Jansen, D.N. Petcher, J.C. Vink, Phys. Rev. D49 (1994) 1606. 

[4] V. Furman, Y. Shamir, Nucl. Phys. B439 (1995) 54. 

[5] J. Schwinger, Phys. Rev. 128 (1962) 2425. 

[6] S. Coleman, R. Jackew, L. Susskind, Ann. Phys. 93 (1975) 267; 
S. Coleman, Ann. Phys. 101 (1976) 239. 

7 



[7] K. Jansen, Phys. Lett. B288 (1992) 348. 

[8] M. Creutz, I. Horvath, Phys. Rev. D50 (1994) 2297. 

[9] Y. Shamir, Nucl. Phys. B417 (1994) 167. 
[10] S. Duane, A.D. Kennedy B.J. Pendleton, D.Roweth, Phys. Lett. B195 (1987) 216. 
[11] H. Dilger, PhD thesis, 1993, University Hamburg, preprint DESY 93-181. 
[12] A. Jaster, PhD thesis, 1995, University Hamburg. 



8 



